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We investigate the effects of suppression of single-particle dispersion near the Fermi level on 
the spin gap and the singlet-pairing correlation by using the exact diagonalization method for 
finite-size systems. We consider strongly correlated one-dimensional models, which have flat 
band dispersions near wave number k = it/2, if the interactions are switched off. Our results for 
strongly correlated models show that the spin gap region expands as the single-particle dispersion 
becomes flatter. The region where the singlet-pairing correlation is the most dominant also 
expands in models with flatter band dispersions. Based on our numerical results, we propose a 
pairing mechanism induced by the flat-band dispersion. 
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One of the characteristic features of hole-doped high- 
T c cuprates is the emergence of flat single-particle disper- 
sion near the X points [i.e. k = (tt,0), (O,7t)].0&§ This 
unusual flatness of the dispersiofl.-|relation is ascribed 
to strong electronic correlations.oB' A recent numeri- 
cal study on the two-dimensional (2D) Hubbard model 
showed that the flat single-particle dispersion near the 
X points is. related to the universality class of the Mott 
transition.© Pairing instabilities ofpSMStems under this 
universality class were also studied, El where the dy- 
namical exponent z changes from four to two for the 2D 
Hubbard and t-J models when a two-particle hopping 
term is switched on. This shows potential instability of 
flat single-particle dispersions towards singlet-pair for- 
mation, if two-particle hopping processes are introduced. 

Taking this into account, we consider a possible sce- 
nario for enhancing the pairing instability: In strongly 
correlated systems, the two-particle hopping process can 
mainly be controlled by the incoherent part of charge 
excitations independently of the coherent part of single- 
particle excitations, while such a coherent single-particle 
process may rather suppress the pairing susceptibility 
through damping and pair breaking processes. There- 
fore, if only the single-particle dispersion near the Fermi 
level is suppressed, the two-particle hopping process de- 
termined from the dispersive band far from_the Fermi 
level would enhance the pairing instabilities.^!! 

In this letter, in order to confirm the occurrence of 
this scenario in a simplified situation, we introduce one- 
dimensional (ID) models ptliat have flat single-particle 
dispersion near half filling.H We clarify how the singlet- 
pairing correlation and the spin gap behave in the pa- 
rameter space of our model. 

We consider the Hubbard model with long-range hop- 
ping terms defined by the following Hamiltonian: 



W H ub = -t^{4a C i+^+ C ^l a C ^) + U ^ n ^ n ii^ ( lb ) 

i,cr i 

H h ± = -t 3 ^(ct CT Ci +3CT + 4+3 CT C ^) 

_ *7 J2( C ^ Cl+7 ° + c \+1o c ^)- ( lc ) 

We take hopping amplitude i, (i = 3, 5, 7) of the long- 
range hopping terms as (i) t 3 /t — 1/3, t 5 /t — tr/t = 0, 
(ii) t 3 /t = 0.5, t 5 /t = 0.1, t 7 /t = 0, and (hi) t 3 /t = 
0.6, ts/t — 0.2, tj/t = 1/35 in order to realize flat disper- 
sion near half filling in the noninteracting case. With this 
tuning, the noninteracting dispersion around go = ±7r/2 
is suppressed up to (q — qo) u , where v = 2, 4 and 6 for 
(i), (ii) and (iii), respectively. The dispersion relations 
of these noninteracting models are shown in Fig. 1. The 
region of flat-band dispersion near half filling is allowed 
to expand further as the hopping range becomes longer. 
Taking the limit of strong correlation (i.e. U — > 00), we 
obtain the following effective Hamiltonian: 



7~i = Ti-t-J + 7Y3_ s itc + Hlrh + HlongJ, (2a) 
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Here, Wi on gj denotes the spin exchange and three-site 
terms originating from long-range hopping. In this letter, 
we mainly study the model defined by the Hamiltonian 
Htm = Tit-. j + + Wi r h for simplicity. Later, we 

discuss the effects of H\ on gj briefly. The operator c^. 
creates an electron with spin a at site i in the subspace 
without double occupancy. The operators Si and rij are 
the spin and number operators at site i, respectively. 
Hereafter, the chain length and the number of electrons 
are denoted by L and N c , respectively. We set t = 1 as 
the energy unit. 

We call models defined by Htjsi with (i) < 3 = 1/3, t 5 = 
t 7 = 0, (ii) tg = 0.5, ts = 0.1, tr = 0, (iii) t 3 = 0.6, t 5 = 
0.2, t 7 = 1/35, and (iv) t 3 = t 5 = t 7 = as tJ3t3, tJ3t5, 
tJ3t7 and tJ3, respectively. The tJ3 model is nothing 
but the effective model of the strong-coiipling limit of 
the Hubbard model up to the order t 2 /Ux3 The ground- 
state properties of this model were imzEstigated in detail 
by Ammon, Troyer and Tsunetsugu.liiP 

The aim of this study is to demonstrate that the sup- 
pression of single-particle hopping processes enhances 
singlet pair correlations. We show this enhancement 
by studying the effects of the above long-range hopping 
terms Wi r h on the spin gap and the singlet-pairing cor- 
relation. 

First, we show numerical results on the phase 
bound ar,y .between the Tomonaga-Luttinger (TL) liquid 
phaseliSEJf and the spin-gap (SG) pliasa. _Wp .emnlny-the 
singlet-triplet level crossing methoo£Zl'll3'oE3'Eil'E3'E3 ) in 
order to determine the phase boundary. The effective- 
ness of this method has been demonstrated through the 
applications. to. ID spin systemsO-HllB' and electron sys- 
tems.E2ltirL3'E3) Using this method—sre determine the 
critical point with high precision,EZrEB' since the finite- 
size correction is .only due to irrelevant fields with scaling 
dimension x — This method is based on the assump- 
tion that the universality class of the phase transition is 
described by th£_level k = 1 SU(2) Wess-Zumino-Witten 
(WZW) modelES 1 In order to check the applicability of 
this method to our models, we have calculated the scaling 
dimension of the lowest singlet excitation (x s ) and that 
of triplet excitation (xt) at wave number k — 2kp. Ac- 
cording to the renormalization-group analysis, the lead- 



ing finite-size corre 
point are written asl 



of x s and Xt near the critical 
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where yo is a coupling constant. As shown in Fig. 2, 
the value of x T = (x s + 3x t )/4 is very close to 1/2 and 
independent of the value of parameter J. This indicates 
that the low-energy physics of our models is described 
as that of the Tomonaga-Luttinger liquid, and that the 
transition belongs to the universality class of the k = 1 



SU(2) WZW model. In fact, as shown in Fig. 3, the 
phase boundaries in clusters of length L = 8, 10, 12, 14 
and 16 are scaled in single lines. This suggests that the 
size dependence is very small. 

Figure 3 shows the phase boundaries between the TL 
liquid phase and the SG phase for the tJ3, tJ3t3, tJ3t5 
and tJ3t7 models. Clearly, the region of the SG phase 
expands further near half filling as the hopping range 
becomes longer. In other words, the suppression of the 
Fermi velocity enlarges the region of the spin-gap phase. 
This may be explained as follows: Generally, low-energy 
properties of fcrmionic systems are mainly determined by 
scattering processes near the Fermi level. The suppres- 
sion of the Fermi velocity would have similar effects as 
the reduction of hopping amplitude t on the low-energy 
physics, because the Fermi velocity is roughly propor- 
tional to t. The long-range hopping terms, which sup- 
press the Fermi velocity, effectively decrease t, while they 
do not alter J and the three-site terms. Then, they in- 
crease the effective J/t for the tJ3 model. In the large 
J ft regime for the tJ3 model, the spin excitation has 
an energy gap, because of the large two-particle hop- 
ping processes such as spin exchange and three-site hop- 
ping.E-f Hence, the suppression of single-particle hopping 
processes due to the long-range hopping terms results in 
the enlargement of the spin-gap phase. 

It should be noted that the spin-gap region extends 
only in the finite-doping regime, not at half filling. This 
is analogous to the behavior of high-T c cunrates and dif- 
ferent iom_that of the dimerized modeLn£ L-thf; t-J-J' 
modelBBEJ and the t-J ladder model |§@0 where 
the spin gap opens not only in the finite-doping regime, 
but also at half filling. This difference results in differ- 
ent classes of proximity effects from Mott insulators and 
hence different pairing mechanisms. Our study provides 
some insight on this issue. 

Next, we consider correlation functions. Generally, in 
ID interacting systems, the long-range behavior of corre- 
lation functions is described by a single correlation expo- 
nent Ifp.EJaEJ The correlation exponent K p is obtained 
as 



K p =7T 



<jQn 2 K 



(4) 



where Drude weight of the ac conductivity <jq and com- 
pressibUity-/t in finite-size systems are respectively de- 
finedBEiy) as 
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A^ c 2 E Q (L,N c + 2) + E (L,N c - 2) - 2E (W 

Here, Eo(<f>) denotes the ground-state energy of the sys- 
tem with twisted boundary conditions with phase factor 
<p. Eq(L, N e ) denotes the ground-state energy of the N e - 
electron system of length L, and n is the electron density 
defined by n = N c /L. The long-range behavior of the 
singlet-pairing correlation function P{r) is expressed in 
terms of K p as P(r) oc l/r 1+1 / Kp if the spin excitation 
is gapless, while P{r) oc l/r 1 / Kp in the spin gap phase. 
If K p is larger than one, the singlet-pairing correlation 
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is the most dominant among other correlation functions. 
We show the contour lines of K p = 1 in Fig. 4. The 
region of K p > 1 expands as the hopping range be- 
comes longer. This means that the suppression of single- 
particle hopping enhances the singlet-pairing correlation. 

Here, we briefhz. refer to the effects of the neglected 
terms in Tii ong j&3 We consider the system defined by 
the following Hamiltonian: 

H = H t -j + W3- S ite + Ti-irh + Ti-iexj , (6a) 

T~L\cxj = -h y^(Sj • Si + 3 — -nin i+3 ) 

i 

+J5^2{Si ■ S i+5 - -nin i+5 ) 

i 

+ <h ^{Si ■ Sj+7 - ^niiii+j), (6b) 

i 

where coupling constant Jj is defined as Jj = J x tj 
(j = 3,5 and 7). The spin-gap phase boundary and 
the contour line of K p = 1 for the above model with 
t 3 = 0.6, t 5 = 0.2, t 7 = 1/35 are shown in Fig. 5. This 
shows that the long-range spin exchange interaction also 
enlarges the spin-gap region and the region of K p > 1. 
The region of singlet pairs may be determined mainly by 
the following two factors: One is two-particle hopping 
processes, which enhance the coherent motion of electron 
pairs and prevent phase separation, and the other is an- 
tiferromagnetic fluctuations, which induce attractive in- 
teractions between up and down spins. For the t-J model 
in the large J regime, the spin-gap phase exists, but the 
antiferromagnetic fluctuations are so large that the spin- 
gap phase is confined in a narrow parameter regime due 
to instability resulting from phase separation.O'EiP For 
the t J3 model, single- and two-particle hopping processes 
are strong enough to prevent phase separationli-P How- 
ever, since the larger single-particle hopping is accompa- 
nied by the three-site term, it requires a large J to reach 
the region of K p > 1. The numerical results shown in 
Fig. 5 indicate that the suppressed single-particle coher- 
ence together with preserved two-particle hopping and 
finite-range spin exchanges contribute to enhancing both 
the effective attractive interactions between up and down 
spins and the coherent pair motion. 

Here, we consider possible realizations that may be 
described by the above models. An example is a three- 
chain system: The central chain is composed of atoms 
with ci-orbitals where the Fermi level is located. The 
other two chains are composed of atoms with 7r-orbitals. 
These 7r-orbitals are designed to lie near the Fermi level 
and to overlap with the d-orbitals of the central chain. 
In such a system, electrons on <i-orbitals can hop to a 
distant site by virtual hopping processes through the 
side chains. The low-energy physics may be described 
by the tJ3t3 model if the parameter is tuned. Another 
is a spiral chain: The orbitals of third neighbor atoms in 
a chain may overlap, if the chain has a spiral structure. 
The overlap between orbitals may produce single-particle 
hopping processes between distant sites. In such a sys- 
tem, the low-energy physics may be described by the 



tJ3t3 model. 

In summary, we have demonstrated that the suppres- 
sion of the single-particle hopping processes near the 
Fermi level enlarges both the spin-gap region and the 
region for the dominance of the pairing. Based on the 
present results, we propose a possible scenario for a pair- 
ing mechanism: If the suppression of the single-particle 
hopping is sufficient, pairing correlations would be en- 
hanced, which may lead to superconductivity. 

Although the models investigated by us are specific 
and somewhat artificial, the underlying physics would 
be universal. The concept of flat-band induced pair- 
ing instability may be appropriate not only for one- 
dimensional models with long-range hopping terms, but 
also for models with short-range hopping terms in higher 
dimensions, if the single-particle hopping processes are 
sufficiently suppressed near the Fermi level close to the 
Mott insulator. 
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Fig. 1. Dispersion relations of noninteracting models with long- 
range hopping terms, where ts/t, t$/t, tj/t take values of (i), (ii), 
(iii) in the text. The horizontal line corresponds to the Fermi 
level at half-filling. For comparison, the dispersion relation of 
the noninteracting model without long-range hopping terms is 
also plotted. 
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Fig. 2. Scaling dimension of the lowest singlet excitation (x B ) and 
that of triplet excitation (xt) at k = 2fcp for the tJ3t7 model at 
quarter-filling in a 16-site cluster. We also plot the quantity 
x- r = (x a + 3x t )/4. 
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Fig. 5. Spin-gap phase boundary and the contour line of K p = 1 
for the tJ3t7 model with the long-range spin exchange term. 
The spin-gap phase boundary and the contour line of K p = 1 
for this model are represented by solid circles and a dash-dotted 
line, respectively. For comparison, we also plot those of the t J3 
model. The spin-gap phase boundary and the contour line of 
K p = 1 for the tJ3 model are represented by open circles and a 
dotted line, respectively. 
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Fig. 4. Contour lines of K p = 1 for the tJ3, tJ3t3, tJ3t5 and 
tJ3t7 models in a 16-site cluster. 



